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ABSTRACT 

Aims. We investigate the cosmological evolution of large- and small-scale magnetic fields in galaxies in the light of present models of formation 
and evolution of galaxies. 

Methods. We use the dynamo theory to derive the timescales of amplification and ordering of magnetic fields in disk and puffy galaxies. 
, Turbulence in protogalactic halos generated by thermal virialization can drive an efficient turbulent dynamo. Results from simulations of 
• hierarchical structure formation cosmology provide a tool to develop an evolutionary model of regular magnetic fields coupled with galaxy 

formation and evolution. 

Results. The turbulent (small-scale) dynamo was able to amplify a weak seed magnetic field in halos of protogalaxies to a few fiG strength 
■ within a few 10** yr. This turbulent field served as a seed to the mean-field (large-scale) dynamo. Galaxies similar to the Milky Way formed 
I their disks at z a 10 and regular fields of yuG strength and a few kpc coherence length were generated within 2 Gyr (at z ~ 3), but field-ordering 
on the coherence scale of the galaxy size required an additional 6 Gyr (at z ~ 0.5). Giant galaxies formed their disks at z a; 10, allowing more 
efficient dynamo generation of strong regular fields (with kpc coherence length) already at z s; 4. However, the age of the Universe is short 
for fully coherent fields in giant galaxies larger than 15 kpc to have been achieved. Dwarf galaxies should have hosted fully coherent fields at 

' Conclusions. This evolutionary scenario can be tested by measurements of polarized synchrotron emission and Faraday rotation with the 
planned Square Kilometre Array (SKA). We predict an anticorrelation between galaxy size and ratio between ordering scale and galaxy size. 
Weak regular fields (small Faraday rotation) in galaxies at z ^ 3 are signatures of major mergers. Undisturbed dwarf galaxies should host fully 
coherent fields, giving rise to strong Faraday rotation signals. 

Key words. Techniques: polarimetric - galaxies: evolution - galaxies: high-redshift - galaxies: interactions ~ galaxies: magnetic fields - radio 
continuum: galaxies 



1. Introduction 

Polarized synchrotron emission and Faraday rotation infeiTed 
the presence of regular large-scale magnetic fields with spi- 
ral patterns in the disks of nearby spiral galaxies (Beck 
2005), which were successfully reproduced by mean-field dy- 
namo theory (Beck et al. 1996; Rudiger & Hollerbach 2004; 
Shukurov 2005). It is therefore natural to apply dynamo theory 
also in predicting the generation of magnetic fields in young 
galaxies at high redshifts. However, there is a lack of obser- 
vational information with which to test models. The available 
data (Kronberg et al. 2008; Bernet et al. 2008; Wolfe et al. 
2008; see also the review by Kronberg 1994), although lim- 
ited, imply that the environments of galaxies are significantly 
magnetized at high redshifts (z ^ 3), with regular magnetic- 
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field strengths comparable to or higher than those at the current 
epoch. Another indication of strong magnetic fields in young 
galaxies is the tight radio - far-infrared correlation, which is 
valid in galaxies at least to distances of z ^ 3 (Seymour et 
al. 2008). To explain the almost linear correlation, the total 
magnetic-field strength must be related to the star-formation 
rate (Lisenfeld et al. 1996, Niklas & Beck 1997). Since obser- 
vations at high redshifts are biased towards starburst galaxies, 
strong magnetic fields presumably do exist in distant objects. 

We now have sufficient evidence that strong magnetic fields 
were present in the early Universe and that synchrotron emis- 
sion should be detected with future radio telescopes such as the 
Square Kilometre Array (SKA). The SKA will spectacularly 
increase the sensitivity and angular resolution of radio obser- 
vations and allow us to observe an enormous number of dis- 
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tant galaxies at similar resolution to that achievable for nearby 
galaxies today (van der Hulst et al. 2004). 

The dynamo was introduced in to astrophysics as a mech- 
anism for transforming the kinetic energy of the motions of 
electrically conductive media into magnetic energy. This con- 
cept describes the amplification of small-scale as well as large- 
scale, regular fields. For practical, mathematical reasons, the 
dynamo equation was separated into a small-scale and a large- 
scale part, two aspects of the same mechanism with widely 
different timescales, but often considered as two independent 
dynamo models. 

The amplification of large-scale magnetic fields and the 
general growth of magnetic energy are not the only merits of 
galactic dynamos. For the evolution of galactic magnetic fields, 
the transformation of a turbulent into a regular field is of sim- 
ilar value. It is the large-scale regular field whose observation 
is most spectacular and whose generation has attracted most 
attention in astrophysics. As far as we know today, the genera- 
tion of regular magnetic fields is a delicate process that requires 
large-scale rotation of the body under consideration. This rota- 
tion produces the a-effect, which amplifies the magnetic field, 
and differential rotation ensures that this amplification is effec- 
tive. This defines the classical "mean-field" galactic dynamo 
as discussed below, taking into account recent results from nu- 
merical simulations of galaxy formation. 

Developing a consistent dynamo model to describe the en- 
tire lifetime of a galaxy requires information about the basic 
physical parameters controlling dynamo action and their evolu- 
tion. Regular large-scale magnetic fields can be generated and 
amplified by the mean-field galactic dynamo in high-redshift 
galaxies, provided that a gaseous, rotating disk already exists. 
The formation of disk galaxies and the epoch of this formation 
are fundamental problems in astronomy. High resolution nu- 
merical simulations of disk formation in galaxies demonstrated 
that a dynamical disk could be formed at redshifts z ~ 5 - 6 and 
even higher (Governato et al. 2007; Mayer et al. 2008). A more 
robust understanding of the history of magnetism in young 
galaxies may help to solve fundamental cosmological questions 
about the formation and evolution of galaxies (Gaensler et al. 
2004). 

In the absence of large-scale rotation, turbulent motions can 
amplify magnetic fields, which are, in general, irregular, small- 
scale fields. This process is therefore known as the turbulent or 
small-scale dynamo (e.g. Zeldovich et al. 1990). The growth 
in magnetic fluctuations is clearly visible in many numerical 
MHD simulations of the interstellar medium (Meneguzzi et al. 
1981; Hanasz et al. 2004; de Avillez & Breitschwerdt 2005; 
Kowal et al. 2006; Iskakov et al. 2007; Gressel et al. 2008; 
Wang & Abel 2008). We note that all numerical simulations 
of the generation of regular magnetic fields in disk galaxies 
begin at the disk formation epoch; we assume that seed fields 
of strength between a few jjG and 10"^ jiG exists already in the 
disk and do not consider the field amplification up to this level. 

High-resolution simulations of protogalactic clouds 
demonstrated that significant turbulence can be generated prior 
to disk formation during the thermal virialization of the halo 
(Wise & Abel 2007). Strong turbulence at these early epochs 
may drive a turbulent dynamo, which can effectively amplify 



the halo magnetic fields. The amplified magnetic fields can 
form the basis of a seed field for the generation of regular 
fields in the disk of a newly formed galaxy. 

The first model of dynamo action in young galaxies by 
Beck et al. (1994) proposed a two-stage dynamo, a turbulent 
dynamo to amplify a weak, seed field in a disk within less than 
a Gyr and a large-scale dynamo to amplify the regular field 
within a few Gyr However, some of the assumptions, e.g. the 
epoch of disk formation, the large disk thickness and the con- 
stant radius, are in conflict with the modern scenario of galaxy 
formation. 

Dynamo theory has been well developed as a branch of 
MHD and, as a first step, we exploit mainly conventional ap- 
proaches to galactic dynamos (see e.g. Ruzmaikin et al. (1988), 
Beck et al. (1996) and Rudiger & HoUerbach (2004)) rather 
than more recent findings in dynamo theory (see e.g. Tobias & 
Cattaneo 2008). In this paper, we develop a simplified model 
for the evolution of magnetic fields in both protogalactic halos 
and galaxies based on the recent numerical developments in 
the study of the formation and evolution of galaxies during the 
epoch of hierarchical structure formation. We use the known 
timescales of the turbulent dynamo and mean-field dynamo in 
puffy and disk galaxies to explore the earliest generation and 
evolution of regular magnetic fields in halos and to describe 
the main phases of magnetic-field evolution in dwarf and disk 
galaxies, as well as the influence of merging events on the mag- 
netic evolution. In Sect.|2l we derive the timescales of field am- 
plification and field ordering in disk and spherical galaxies by 
the mean-field dynamo. The timescale of the turbulent dynamo 
is discussed in Sect.[3l and the origin of seed magnetic fields in 
Sect, m In Sect. |5] we summarize developments in the forma- 
tion and evolution of disk and dwarf galaxies, and in Sect. |6l 
we develop an evolutionary model of regular magnetic fields 
in galaxies. Perspectives on detecting polarized radio emission 
from nearby and from distant galaxies are discussed in Sect. |2l 
and conclusions are drawn in Sect. [8] 

Throughout the paper, a flat cosmology model (Qa + Q.m - 
1) is used with Q.„, - 0.3 and Hq - 70 kms"' Mpc"'. 

2. Magnetic-field generation in galaxies by the 
mean-field dynamo 

In this Section, we refer to thin-disk galaxies that rotate dif- 
ferentially. Theoretical analysis of dilTerentially rotating disks 
demonstrated that a disk is only locally stable to axisymmet- 
ric perturbations if the Toomre's stability criterion is fulfilled 
(Julian & Toomre 1966), i.e. only cool disks in which both the 
gas and stars have low velocity dispersion can survive for a 
long time (Binney & Tremaine 1987). Hence, the velocity dis- 
persion in disk galaxies at higher redshifts is probably similar 
to that in nearby galaxies. In the following, we adopt the usual 
values of v = 10 km s"' and a scale of turbulence of / - 100 
pc. Turbulence in galaxy disks is assumed to be driven by su- 
pernovae (Ruzmaikin et al. 1988; Korpi et al. 1999; de Avillez 
& Breitschwerdt 2005; Gressel et al. 2008). The other basic 
parameters of the galactic mean-field dynamo are the angular 
velocity of the galactic rotation Q, the disk thickness h, the disk 
radius R, and the gas density p. 
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2.1. Timescales of regular magnetic fields in disk 
galaxies 

The conventional theory of the galactic mean-field dynamo in 
disk galaxies allows the generation of a regular galactic mag- 
netic field as a result of joint action of diff'erential rotation 
D. and helical turbulent motions of interstellar gas. The lat- 
ter is responsible for the so-called a-eff'ect (a). According to 
Ruzmaikin et al. (1988), the intensity of both generators can be 
characterized by two dimensionless numbers, and Ra'. 



R,,, — 



1/ h 



/?\500pc/ \20kms-ikpc" 
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where Q and h are the angular velocity of the galaxy and the 
half-thickness of a galactic disk (normalized to the galactic an- 
gular velocity and half-thickness, respectively) and fi is the tur- 
bulent difFusivity. A flat rotation curve (rdQ./dr ~ Q) was as- 
sumed in deriving the above relations. The parameters a and p 
cannot be measured directly from observations, and our aim is 
to reduce them to observable parameters. According to simple 
mixing length theory the turbulent difflisivity is 



Iv 



(3) 



where I is the basic length scale of turbulence, v is the turbu- 
lent velocity, and the numerical factor characterizes the dimen- 
sions of the space. Typical numbers, as found from observa- 
tions of nearby galaxies, are / - 100 pc and v = 10 km s"' 
(e.g. Ruzmaikin et al. 1988). The a-efl'ect is determined by the 
Coriolis force and density gradient (Krause & Radler 1983): 



a = 



(4) 



The joint action of both generators can be described by a 
so-called dynamo number: 



m^Ra>Ra^9i — \ 



(5) 



The absolute value is taken because is usually considered to 
be negative. The dynamo growth rate (F) is given by 



(6) 



where is the time of turbulent diffusion. The dynamo 
timescale for amplification of the regular field is /* = F ' (e- 
folding time), and we derive 



^ ~ Til 



(7) 



from Eqs. Q, (|5]l and (|6|l. This estimate is the result of com- 
bining standard estimates from galactic dynamo theory, which 
has not been presented before. 



During its lifetime T, a galaxy completes rotations. 
Assuming Q » 20kms"' kpc"', T - 10'" years and hjl x 5, 
we derive a; 30 and t* * STjilnN) ~ 2.5 x 10^ yr. Hence, the 
conventional galactic mean-field dynamo can amplify the reg- 
ular magnetic field by a factor of at most 10'^, while in practice 
it can be lower (see Sect. 4). 

The galactic dynamo is a threshold phenomenon. The reg- 
ular magnetic field is only maintained (or amplified) if the 
dynamo number is larger than a critical value |Dd| > |£>cr,dl- 
Ruzmaikin et al. (1998) derived |Z)cr,dl - 7 from the numeri- 
cal simulations of galactic dynamo models, which they deter- 
mined to be the minimal value required for the dynamo effect 
to overcome turbulent diff'usion. If |Z)d| < |Z)cr,dL the regular 
magnetic-field decays. The typical value of dynamo number 
for the Milky Way (called MW hereafter) is |Dd| - 9 (note that 
/jQ/v ~ 1 for the MW), indicating that the dynamo mechanism 
acts in the MW. If /lO/v I, a disk galaxy cannot amplify 
or maintain regular magnetic fields by means of the mean-field 
dynamo mechanism. 

The expression for t* (Eq. |7]l can be written in terms of the 
magnetic-field strength. The time required to amplify the large- 
scale magnetic field from Bq to B, is given by 



f = f. ln(B./Bo); 



(8) 



lfi>T, the dynamo is in its kinematic regime until the present 
time. Hence, the minimum seed field required for dynamo ac- 
tion is 



Bo = B*exp(-r/f*). 



(9) 



We note that f is the e-folding timescale for amplification 
of the magnetic field. The resulting magnetic field is spatially 
ordered on a coherence scale Ic of a few kiloparsecs (kpc), de- 
pending on the disk half-thickness. Full ordering of the field on 
the coherence scale of the radius R (/c = /?) of a galaxy takes a 
far longer time. According to Moss et al. (1998), this ordering 
timescale t is 



t - 



Vl^ P h 



-1/4 



(10) 



Magnetic-field ordering on the scale of /c is not sim- 
ply an instability that can be adequately described by an e- 
folding time but a more delicate process called "front propa- 
gation" (the so-called Kolmogorov-Petrovsky-Piskunov effect, 
Kolmogorov et al. 1937; for review in the context of dynamo 
theory, see also Zeldovich et al. (1990)). For a crude estimate, 
one can use t from Eq. ( fTOl i in the sense that, for times t < t, 
the ordering remains in progress, while for t > t, the ordering 
process is already finalized. 

Equation ( fTOb can be presented in terms of diffusion 
timescales in the vertical (T/,) and radial (Jr) directions. The 
timescale for full ordering t{lc - R) is far longer than the am- 
plification timescale f (see factor R/h in Eq.llOll as well as the 
diffusion timescale in the vertical direction: 



rA = ^-7108yr, 



(11) 



4 



T. Arshakian et al.: Evolution of magnetic fields in galaxies 



and t <K Tr, 



= — ^ 400 Th » T. 



(12) 



The ordering timescale is governed by the geometrical mean of 
both timescales 



(13) 



For MW-type galaxies, t is only a few times shorter than 
the galaxy lifetime T, indicating that transient field configura- 
tions, substantially different from the leading eigenmode of the 
galactic dynamo, can survive in some galaxies until now. This 
may explain e.g. the large-scale field reversals observed in the 
MW. 

After the ordering timescale f, a steady-state magnetic field 
is maintained. The configuration of this field is expected to be 
similar to the leading eigenmode of the galactic dynamo, satu- 
rated by nonlinear effects (see Sect. l2.2l l. 

2.2. Equilibrium magnetic-field strengtli in disk 
galaxies 

The initial growth of the regular magnetic field due to the galac- 
tic dynamo is saturated at some level, which is known as the 
equilibrium magnetic-field strength (Beq) at which dynamo ac- 
tion reaches equilibrium with turbulent dissipation. The equi- 
librium magnetic-field strength is related to the equipartition 
field strength (B,): 



2 
or 

B, -. 



Bl 



4np. 



(14) 



(15) 



The more elaborated estimate by Shukurov (2004) takes 
into account the intensity of dynamo action measured by Dj: 



Beq — B 



Da 



- 1 . 



Using Eq. ([Sj and D„ - 7, we obtain 



(16) 



eq 



4np 1.281 — 



1 , 



(17) 



The concept that Beq is determined by B» is an oversimpli- 
fication. A deeper understanding of dynamo saturation based 
on magnetic helicity conservation (e.g. Kleeorin et al. 2002) 
yields a similar estimate for Beq- 

We conclude that the timescale of dynamo growth and equi- 
librium field strength are determined by different governing 
parameters. In particular, if v becomes larger (and the other 
parameters are fixed), the dynamo timescale remains constant 
(Eq.|7J, while the equilibrium strength increases. On the other 
hand, the dynamo number strongly decreases with increasing v 
(Eq.|5]l and may drop below the critical value. 



2.3. Timescale of regular magnetic fields in 
quasi-spherical galaxies 

We estimate the dynamo timescales in spherical or "puffy" 
galaxies with li/R > 0.1. For a spheroid, h ^ R: 



m = RaR^ - 9 



(v 



while the dynamo growth rate is 



(18) 



(19) 



(Sokoloff 2002). This yields the following estimate of the dy- 
namo timescale t** of a spheroid: 



_ 3 ( V y/3 R 



9^-'^\RD.I ill' 



(20) 



We note that the final term is the dynamo timescale for a disk 
f*, where h is replaced by R. Assuming that v = 10 km s"' and 
RQ. - 250 km s we estimate that 



f" = 0.23t*(R), 



(21) 



where f(R) means that in Eq. ^ his replaced by R. Since R is 
generally a few times larger than h, the timescales for disk and 
spherical galaxies are more or less comparable. 

The dynamo number of a quasi-spherical galaxy is calcu- 
lated from the ratio of Eq. to Eq. (HI, = \Di\(R/hf = 
10^ - 10^ for R/h = 10 and |Dd| = 10 - 100 (e.g. Belvedere 
et al., 1998). The critical dynamo number |Z)cr,sl for a quasi- 
spherical body is within the wide range between 300 and 5x10^ 
(cf. Sokoloff & Shukurov (1990) and Sokoloff et al. (2008)), 
depending on the details of the rotation curve and the spa- 
tial distribution of the dynamo-governing parameters. For the 
above estimates, IDJ ~ 5000. In this paper, we assume that 
\Df,\ > |Z)ci,sl and do not specify the suppression of the mean- 
field dynamo in spherical bodies. 

3. Magnetic-field generation by the turbulent 
dynamo 

Magnetic-field generation by the turbulent (small-scale) dy- 
namo in galaxies requires neither large-scale rotation nor a 
disk, only turbulence. It produces magnetic fields on scales 
comparable to the basic scale of galactic turbulence on 
timescales far shorter than that for the conventional mean-field 
galactic dynamo. The timescale (e-folding time) of the small- 
scale dynamo is given by 



(22) 



(see Batchelor (1950), Landau & Lifshitz (1959), and 
Kazantsev (1967)), while the equilibrium strength b remains 
comparable to B,. 

Considering that interstellar turbulence produces random 
small-scale magnetic fields as well as regular large-scale ones, 
we compare their strengths. We note that the term -\/D/Dci - 1 
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in Eq. (fTST i is usually less than unity, so that the ratio of small- 
scale and large-scale field strengths /j/Beq is expected to be 
slightly larger than unity. A ratio of ^j/Beq - 1-7 is estimated 
from observations (Ruzmaikin et al. (1988), see however Beck 
et al. (2003) for a systematic bias imposed on b/B^q). For the 
MW, this ratio, derived from the above theoretical arguments, 
is 

^ - , ^ - 1.9, (23) 

^eq ^1,28(^)2-1 

in good agreement with the observational value. 

4. Origin of seed magnetic fields 

The small-scale and large-scale dynamo mechanisms cannot 
explain the origin of magnetic fields. Seed magnetic fields are 
required for the dynamo mechanism to amplify and maintain 
the magnetic field themselves. One possible origin of a seed 
field is the Biermann battery mechanism, which was initially 
suggested in a cosmological context by HaiTison (1970) and ap- 
plied to protogalaxies by Mishustin & Ruzmaikin (1971) (see 
Ruzmaikin et al. (1988) for a review). The typical magnetic- 
field strength produced by the battery mechanism is estimated 
to be between Bseed = lO^^o ^^^j jg^^a q 

In the context of a hierarchical structure formation scenario, 
the Biermann battery is able to generate a field of order 10"* G 
in a protogalactic cloud at z ^ 40 (Pudritz & Silk 1989; Davies 
& Widrow 2000). This theoretical prediction is consistent with 
magnetohydrodynamical simulations of Population 111 star for- 
mation by Xu et al. (2008). Magnetic fields of order ^ 10"'^ 
G are generated predominantly by the Biermann effect, early 
in the evolution of the halo (z - 40) at low densities and large 
spatial scales (400 pc). 

Other possible mechanisms for generating magnetic seed 
fields are kinetic plasma instabilities, the so-called Weibel in- 
stability (Weibel 1959) in cosmological shocks (Schlickeiser & 
Shukla 2003; Medvedev et al. 2004) or cosmological perturba- 
tions (proposed by Takahashi et al. (2006)) that may produce 
field strengths of ~ 10 G (see also the review by Semikoz & 
Sokoloff^(2005)). 

In summary, magnetic fields can be generated by various 
mechanisms and during various stages of cosmological evo- 
lution. It is, however, difficult to generate a regular magnetic 
field on galactic scales that is of reasonable strength in the 
pre-recombination era when the Universe was homogeneous 
and isotropic; the creation of magnetic fields instead requires 
timescales and peculiar motions with respect to the homoge- 
neous and isotropic cosmological background that are incon- 
sistent with the assumptions of isotropy and homogeneity. 

The variety of possible mechanisms for generating seed 
magnetic fields increases with the development of spatial struc- 
tures in the Universe. The first magnetic fields that represented 
seeds of galactic magnetic fields were probably created some- 
time between the epoch of reionization and formation of galax- 
ies (10 ^ z ^ 40). 
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A mechanism creating small-scale magnetic fields with a 
vanishing mean value is able to randomly generate a weak 
large-scale field component of 

= /,A?-i/2 (24) 

where is the number of turbulent cells in a galaxy. This is suf- 
ficient to represent a seed for the large-scale galactic dynamo, 
as pointed out by Ruzmaikin et al. (1988), Beck et al. (1996), 
and Subramanian (1998). The key point is that the number of 
turbulent cells is large but not enormous, much less then the 
Avogadro number A = 6 x 10^^, the typical number of particles 
in laboratory statistical physics, where terms of order N^^^^ can 
be neglected. 

Using = 7TR^2h/(4Kp/3) a; 8 x 10* for a disk galaxy with 
R - 10 kpc, h = 500 pc and / = 100 pc and assuming that the 
small-scale magnetic field is on the equipartition level, b ^ S„ 
we obtain 

B, * B,/300. (25) 

Hence, the large-scale dynamo has to amplify the seed field by 
only a factor of about 300, much less than the upper limit of 
the amplification factor of 10'^ obtained from Eq. (7), and the 
corresponding amplification timescale is then far shorter than 
the galaxy age. 

The evolution of galactic magnetic fields initiated with a 
seed field generated by the Biermann battery (in the epoch of 
protogalaxy formation of protogalaxies), which was followed 
by field amplification by the turbulent small-scale dynamo, and 
further amplification and ordering by the large-scale dynamo. 
A deeper understanding of the process would also require in- 
clusion of the magneto-rotational instability (Kitchatinov & 
Riidiger 2004), although this is beyond the scope of the present 
paper. 

5. Formation and evolution of disk galaxies 

We review results of recent simulations of the hierarchical 
merging of dark halos, and the formation of protogalaxies and 
their evolution (no concise review of the recent numerical re- 
sults exists so far). We attempt to identify the mechanisms re- 
sponsible for generating and amplifying turbulent and regular 
magnetic fields in the halo and disk of galaxies, in developing 
a simple model for the evolution of magnetic fields and deter- 
mining the earliest epochs at which large-scale regular mag- 
netic fields are formed and amplified in disk galaxies, of the 
type considered in the next Section. The formation and evolu- 
tion of regular large-scale magnetic fields is related intimately 
to the formation and evolution of disks in galaxies in terms 
of geometrical and physical parameters such as the radius R, 
half-thickness h, and angular momentum Q. of the galaxy, tur- 
bulence velocity v, turbulence scale I, and density p of the gas. 
These parameters evolve dififerently for different galaxy types 
and sizes. 

We adopt the standard hierarchical cold dark matter cos- 
mology in which the cosmic structures and galaxies assem- 
ble by the merging of small dark matter halos (Kauffmann 
et al. 1993; Baugh et al. 1996). The first very massive stars 
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100 Mq) are predicted to have formed in dark matter halos 
of virial mass ~ 10'' Mq at redshifts z ^ 20 (e.g. Tegmark et al. 
1997; Omukai & Palla 2003). The dark matter halos merged, 
assembling the first protogalaxies of masses ^ 5 x IO^Mq at 
Z ^ 10 (Press & Schechter 1974; Rees & Ostriker 1977). The 
gas and stellar components of protogalaxies may have had a 
quasi-spherical or disk-like distribution, depending on the ini- 
tial conditions of simulations (Brook et al. 2004; Wise & Abel 
2007; Greif et al. 2008). The first galaxies were able to re- 
tain photo-heated gas and maintain self-regulated star forma- 
tion (e.g. Ricotti et al. 2008), and developed into disk galaxies. 

In this scenario, we identify two main cosmological epochs 
at which magnetic fields can be generated and amplified by 
different dynamo mechanisms. In the first epoch of the viril- 
ization and merging of dark matter halos, which occurred be- 
tween 35 ^ z ^ 10, the mass of protogalaxies assembled via 
gas accretion onto the halos and minor/major mergers of dark 
matter halos. High-resolution simulations of early pre-galactic 
halos (Wise & Abel 2007; Greif et al. 2008), which include 
primordial gas cooling and mass accretion history, showed that 
significant turbulence was generated in all cosmological halos 
during thermal virilization. The gas accreted from the IGM was 
shock-heated to the virial temperature. This mode of accretion 
(hot accretion) worked effectively at high redshifts in low-mass 
halos and generated turbulence on scales comparable with the 
size of the infalling gas. As radiative cooling became efficient, 
the gas attempted to virialize by increasing its kinetic (turbu- 
lent) energy, which was achieved by radial infall and turbulent 
motions. The turbulence became supersonic with Mach num- 
bers ranging from 1 to 3 (Wise & Abel 2007). 

In this mode (cold accretion), massive filaments could form 
molecules, which allowed efficient cooling of the filaments and 
their flow at high velocities (~ 20 km s"') towards the center 
of the protogalaxy (Greif et al. 2008). Inflows of cold gas were 
supersonic with Mach numbers ~ 10 and reached the central 
part of the galaxy, generating significant turbulence with Mach 
numbers ranging from 1 to 5. 

Virial turbulence may have been most important in halos 
with efficient radiative cooling rates, for example in halos of 
masses below 10'^ Mq that could be cooled effectively by Lyo- 
emission (Birnboim & Dekel 2003; Dekel & Birnboim 2006). 
For low-mass protogalaxies, the largest driving scale of turbu- 
lence was ~ 200 pc (1/3 of the virial radius) and the r.m.s. 
velocities was around 20 km s"' (John Wise, private commu- 
nication). The dominant role of virial turbulence played an im- 
portant role in star formation in the regions with density en- 
hancements and the amplification of magnetic fields by means 
of the turbulent dynamo (see Sect. 16. lb . 

Mergers could also have generated turbulence that started 
on length scales of merging halos and developed to small scales 
by means of turbulence cascades. Another driver of turbulence 
was the Kelvin-Helmholtz instability between bulk ffows and 
virialized multi-phase gas (Takizawa 2005), which may have 
occurred during minor mergers. Simulations of protogalaxies 
in the redshift range between 35 and 15 demonstrated that the 
turbulent Mach number depends on merger history. The Mach 
number was supersonic, reaching values between 2 to 4 in the 
case of two subsequent major mergers. Before the central col- 



lapse occurred in these halos, turbulence was mostly subsonic 
(< 1 ; Wise et al. 2007). Major mergers of halos resulted in high 
turbulence Mach numbers and hence high turbulent velocities. 

The beginning of the second epoch manifested itself by the 
formation of an extended large-scale disk in the first galax- 
ies. Highly resolved simulations showed that the first galaxies 
were born at redshifts z ^ 10, after the atomic cooling crite- 
rion was fulfilled and they were able to retain the photo-heated 
gas (Greif et al. 2008). Brook et al. (2004) demonstrated that 
the thick stellar disk was created during an epoch of multiple 
mergers of gas-rich halos, which were abundant at high red- 
shifts. The population of stars forming during the merging pe- 
riod formed the thick disk. The angular momentum of a signif- 
icant fraction of the accreted gas (i.e. the angular momentum 
of the protogalactic halo) resulted in the rotation and flattening 
of the formed galaxy. The forming thick gas disk was dynam- 
ically hot, resulting in the formation of high-velocity stars in 
the thick disk. After the epoch of violent merger events, the gas 
was accreted and formed a smooth thin disk. 

Recent models of disk formation have included hierarchical 
dark matter cosmology, gas cooling, star formation, and super- 
nova (SN) feedback (Governato et al. 2007; Stringer & Benson 
2007) and benefitted from the increased resolution of numerical 
simulations (Mayer et al. 2008). High mass resolution and spa- 
tial resolution are important in simulating large MW-type disks 
(Mayer et al. 2008). The formation of disks at high redshifts 
(z > 2) can be simulated with relatively high spatial and mass 
resolutions (Governato et al. 2007). The disk started forming in 
galaxies with M 10" Mq immediately after the last gas-rich 
major merger with substantial stellar feedback (Elmegreen et 
al. 2005; Springel & Hernquist 2005) at ziim ~ 2. Preserving 
the angular momentum from the merger, the gas cooled and 
formed the large-scale exponential disk. This effect was more 
apparent in numerical runs with stronger feedback and higher 
resolution. After the stellar disk had formed, it remained fairly 
stable against merger events provided there was sufficient gas 
to form stars (Springel & Hernquist 2005). As discussed by 
Kaufmann et al. (2007), internal and external UV background 
radiation would prevent cooling below ~ 10"* K in the warm 
galactic gas by providing pressure support dynamically com- 
parable to angular momentum support. They showed that pres- 
sure support was less important in large halos (^ 10^ Mq). 
This created a tendency to form a thin disk in massive, isolated 
galaxies, while their less massive counterparts (dwarf galax- 
ies) remained spheroidal, puffy systems with more gas and less 
efficient star formation than in the larger MW-type galaxies. 
The disk could settle down (and survive until the present time) 
in massive halos after the epoch of gas-rich mergers, even at 
redshifts z S 5 (Lucio Mayer, private communication). Future 
higher-resolution numerical simulations will allow to trace the 
disk formation back to z ~ 10. 

The advanced disk-formation model presented by Stringer 
& Benson (2007) showed that the resulting disks of the iso- 
lated galaxy (no merger) and the galaxy formed through mi- 
nor mergers had very similar evolutionary behaviors. The disk 
component was not destroyed and evolved into a disk galaxy. 
The morphology of the formed disk galaxies at z ^ 1.5 re- 
mained unchanged in the case of no merger or only minor 
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mergers, while in the case of frequent minor mergers or a major 
merger the disk galaxy could evolve into a disk galaxy having 
a spheroidal or bulge component, or into an elliptical galaxy 
without disk (see Sect. l6.4l i. 

If a late major merger destroyed the stellar disk at later 
epochs (z ~ 1), the gaseous and stellar disks could recover, 
if there was sufficient gas remaining to form new stars. 

It is generally believed that, in the hierarchical cold dark 
matter Universe, the massive disk galaxies were formed first, 
while the majority of present-day galaxies along the Hubble se- 
quence such as elliptical galaxies and disk galaxies with bulges 
formed at later epochs by the merging of giant disk galaxies 
(Kauffmann et al. 1993; Baugh et al. 1996). The existence of 
massive galaxies at high redshifts is verified by the detections 
of disk-like or clumpy disk galaxies between 1.4 < z < 3 
from deep near-infrared imaging (Labbe et al. 2003) and at 
z ~ 5.5 (Elmegreenet al. 2007). Observations of luminous star- 
forming galaxies at high angular resolution showed that large 
and massive rotating protodisk galaxies were present already at 
z ~ 2 to 3 (Genzel et al. 2006). 

6. Evolution of magnetic fields in galaxies 

6.1. Three-phase model 

In the hierarchical formation scenario, we identify three main 
phases of magnetic -field evolution in galaxies. In the first 
phase, the seed magnetic fields of order ^ 10"'^ G (Pudritz & 
Silk 1989; Davies & Widrow 2000; Xu et al. 2008) were gener- 
ated in dark matter halos by the Biermann battery mechanism, 
well before the formation of first massive stars at z ~ 20. 

The second phase started at z ~ 20 (or f ~ 0.5 Gyr) when 
the merging of halos and virialization generated turbulence in 
the halo. According to Wise & Abel (2007), this epoch was 
dominated by turbulence generated during the thermal virial- 
ization of halos. There the turbulent dynamo could effectively 
amplify the seed magnetic field of halos to the equilibrium level 
(Eq.l23Tl. Bs ~ 20 ^G, on timescales given by Eq. (l22l i. 

In the third phase, the mean-field dynamo mechanism 
started acting in the newly formed galaxies at z ~ 10. The first 
gas-rich massive galaxies (^ 10^ Mq) formed extended thin 
disks {hjR £ 0.1) at Zdisk - 10 after major merger events, as 
discussed in Sect. |5] A weak, large-scale magnetic-field com- 
ponent was generated in the disk from small-scale magnetic 
fields of the halo (Eq. l24l l. amplified to the equilibrium level 
in the second phase. Then, the "disk" mean-field dynamo am- 
plified and ordered weak, regular magnetic fields on timescales 
given by Eq. (|7]i and Eq. (fTOl i. respectively. If the formed disk 
was thick {hjR ^ 0.1) or the disk had not formed, the "quasi- 
spherical" mean-field dynamo acting in puffy objects would be 
switched on. It would amplify the field to the equilibrium level 
and order the magnetic fields on scales of a few kiloparsecs 
and on timescales given by Eq. (|20] | and Eq. ( fTOb {h x R), re- 
spectively. If a thin disk in puffy objects had been formed at 
later epochs, the "disk" mean-field dynamo would have domi- 
nated and amplified the regular fields at the equilibrium level. 
This scenario is more appropriate for dwarf galaxies 10^ 
Mq) that could not maintain the disk at early epochs because 



of strong UV background radiation, and may have formed a 
disk at later epochs (z ~ 2) when UV radiation was weak and 
allowed efficient radiative cooling. 

If the disk was destroyed by a late gas-rich, major merger 
event, the regular magnetic field must have dissipated and in- 
creased the amplitude of the turbulent magnetic field. During 
the disk recovery, the regular magnetic field must have been re- 
generated and amplified back to its equilibrium level (Eq.ll7b. 
The seed field for the large-scale dynamo is taken from Eq. ( l24l l 
using according to the appropriate galaxy radius. 

The sketch of the evolution in the amplitude and ordering 
scale (Eq.fTOli of magnetic fields in dwarf and spiral galaxies is 
presented in Figs. 1-3. 

6.2. Assumptions of the model 

We assumed an average gas density of a halo (adiabatic model) 
within its virial radius of about 600 pc atz - 15ofs! 6x 
10"^"^ gcm"^(Wise & Abel 2007). The gas densities of nearby 
disk galaxies are calculated to be in the range of between 
10"^"* gcm""^ and 15 x 10"^^ gem"-' for a random magnetic- 
field strength in the range ~ (2 - 10) fiG (Beck et al. 1996). 
For simplicity, we assume ffie same gas density for halo- and 
disk-dominated galaxies and that their gas density and angu- 
lar rotation velocity are unchanged during the evolution, i.e. 
p - 10"^^ gcm"^ and Q = 20kms"' kpc"'. 

Halos. We assume that virial turbulence dominates during 
the epoch of merging and virialization of dark matter halos 
(Wise & Abel 2007). The typical turbulence velocity and the 
largest scale of turbulence, 20kms"' and 200pc in halos of 
masses ^ 10^ M0 (John Wise, private communications) and 
are adopted for halos of masses 

Galaxies. For massive and MW-type disk galaxies (^ 
10'° Mq), we assume that (a) the gaseous and stellar disks 
formed at Zdisk ~ 10; (b) the turbulence scale and turbulence 
velocity were driven by SN explosions and had values close to 
those observed in present-day disk galaxies such as the MW, 
/ = lOOpc and v = lOkms '; (c) the scale height of the disk, 
h = 500 pc, remained unchanged as demonstrated by N-body 
simulations of late-type disk galaxies at z < 1 (Brook et al. 
2006), while its scale length in the plane decreased according to 
R,^R (1 -i-z) ""^^ (Trujillo et al. 2006), where R(z=0) = lOkpc. 
This is in reasonably good agreement with the observed disk 
thickness ratio of h/R - 0.08 at z ~ 1 (Reshetnikov et al. 2003) 
and h/R ~ 0.15 at higher redshifts (Elmegreen et al. 2005). 

For typical low-mass puffy galaxies (^ 10'° Mq), we adopt 
R = 3kpc, h K R ^ 3kpc, and Q = 20km s ' kpc 
The ordering timescale of ffie field in these galaxies is calcu- 
lated from Eq. ( fT3] ) using h ^ R. We assume that the "quasi- 
spherical", mean-field dynamo amplified the regular magnetic 
field in puffy galaxies {h/R > 0.1), while ffie "disk" mean-field 
dynamo was effective in galaxies with a thin disk {h/R < 0.1). 

6.3. Isolated disk galaxies 

We consider the evolution of magnetic fields in giant disk 
galaxies (called GD hereafter; /?(z=0)= 20kpc), MW-type 
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galaxies (lOkpc), and dwarf galaxies (called DW hereafter; 
^(z=0)= 3kpc). We assume that the seeds of turbulent mag- 
netic fields of strength ^ 10 '**G existed in the protogalaxies 
of present disk galaxies, at z » 35. Virial turbulence could am- 
plify turbulent magnetic fields in merging dark-matter haloes 
via the small-scale dynamo from 10 G seed field strength to 
reach the equipartition field strength of 2.2 x 10"^ G at z ^ 11 
(0.4 Gyr) within a short period of time, « 3 x 10'^ yr (Eq. |22l 
see Figs. [T] and |2]l. For simplicity, we consider that all types of 
galaxies reached the equipartition level of turbulent fields at the 
same epoch, z ^ 1 1 . 

MW-like galaxies. A sketch of the evolution in the magnetic 
fields of isolated MW-type galaxies with Mmw ~ 1O"M0 is 
shown in Figs.[T]and|2] At Zdisk ~ 10, the disk was formed and 
evolved in isolation (no major merger) to the present time. The 
size of the galaxy was Rio(z = 10) = R(l+zy'^'^^ - 3.4kpc and 
h/Rio = 0.14 > 0.1, and the "quasi-spherical" mean-field dy- 
namo amplified the field until z « 4, at which time, became 
less then 0. 1 . The "disk" mean-field dynamo then became sig- 
nificantly more important and amplified the regular large-scale 
field within a: 1 .5 Gyr, and reached its equilibrium state at z ~ 3 
(Fig ID- At this epoch, the regular field was ordered on a scale 
of a few kiloparsecs (Fig. [3]) and finally reached a coherence 
scale similar to the size of the MW (lOkpc) at z ~ 0.4. 

The earliest regular magnetic fields of equipartition 
strength were generated in fd.giobai ~ 1 .4 Gyr after disk forma- 
tion (Fig.[T]i, while the ordering of magnetic fields on a length 
scale similar to that of a MW-type galaxy was complete after 
~ 9 Gyr (Fig. |3]l. Hence, present-day, MW-type galaxies are 
expected to host fully ordered regular fields. 

This evolutionary scenario also coiTesponds to the evolu- 
tion of a disk galaxy that experiences multiple minor mergers 
that do not destroy its disk. Multiple and even multiple fre- 
quent minor mergers may have caused local distortion of the 
disk, resulting in a high star-formation rate that shifted the start 
of large-scale dynamo action to a later epoch. This will be dis- 
cussed in more detail in the next Section. 

If amplification by the turbulent dynamo is not consid- 
ered, the mean-field, large-scale dynamo would have to am- 
pUfy a seed field of 10 G, and the time to reach the equilib- 
rium strength of the regular magnetic field would be far longer 
(=s 9 Gyr, until a redshift of z ~ 0.3) (Fig.[T]l. However, this sce- 
nario is unrealistic because sufficient turbulence exists in young 
galaxies to drive the turbulent dynamo. Even without turbu- 
lence, another mechanism may enhance the seed field strength 
by a factor of <: 10'^: the dissipative collapse into a disk (Lesch 
& Chiba 1995) at z ~ 10. This would enable the equiparti- 
tion level to be reached earlier, within » 7.5 Gyr In case of 
turbulence, this process is not required because the turbulent 
field is already strong at this redshift. However, if the collapse 
of the halo occurs at z ^ 10, when the field amplitude remains 
smaller than the equipartition level (see Fig.lTJ, this mechanism 
may further accelerate the field amplification by the turbulent 
dynamo. In the following, we only consider amplification of 
the seed field by the turbulent dynamo. 

Giant galaxies. Late-type disk galaxies of a disk scale 
length S 10 kpc are rare as verified by analyses of I-band 
images (Courteau et al. 2007). In these giant disk galaxies 
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Fig. 1. Evolution of magnetic fields in MW-type disk galax- 
ies: magnetic-field strength versus cosmic epoch. Evolution of 
the small-scale magnetic field generated by the turbulent dy- 
namo (thick dashed-dot-dot-dot line) and the large-scale mag- 
netic field generated by the large-scale mean-field dynamo in 
quasi-spherical galaxies (thick solid line) or in thin -disk galax- 
ies (thick dashed-dot-dashed). Dissipation of the field because 
of a major merger event is presented by a dotted line. The 
lower curve traces the evolution of regular magnetic fields gen- 
erated by the pure large-scale dynamo mechanism (no amplifi- 
cation by the turbulent dynamo). The two horizontal thin solid 
lines represent the equipartition and equilibrium magnetic-field 
strengths (upper and lower lines, respectively). The vertical 
thin dotted lines indicate redshifts from to 10. 
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Fig. 2. Evolution in the magnetic-field strength of dwarf galax- 
ies (DW), MW-type, and giant galaxies (GD). The meanings of 
different fine types are the same as in Fig. [T] 



(Mgd ~ lO'^M©), the ratio hlRxQ{z = 10) = 0.5/6.8 < 0.1 
implies that the "disk" mean-field dynamo has already been 
switched on at z ~ 10 (Fig.|2]l, had amplified the regular mag- 
netic field within only s; 1 Gyr, and reached the equilibrium 
state akeady at z ~ 4. However, as shown in Fig. [3] the regular 
magnetic fields can be ordered only on scales of 15 kpc until 
the present time. Hence, the mean-field dynamo cannot gener- 
ate coherent magnetic fields over the size of the disks of giant 
galaxies {R ^ 15 kpc). 
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Fig. 3. Scale of ordering of regular magnetic fields with cosmo- 
logical epoch. The evolution of the ordering magnetic fields is 
shown for dwarf galaxies (DW; bottom line) and disk galaxies 
(MW and GD; top Une). 



Dwarf galaxies. We assume that, in low-mass halos, the 
turbulent magnetic fields evolve in the same way as in their 
massive counterparts. N-body/smoothed particle hydrodynam- 
ics simulations demonstrated that low-mass galaxies did not at 
first form thin disks because of the ionizing UV background 
radiation, which prevented cooling of the warm galactic gas 
(Kaufmann et al. 2007). They were more spheroidal, puffy 
systems {hjR « 0.3) of rotational velocities ■x, 40 km s"' and 
mass M » 10'" M0. In thick galaxies, the "quasi-spherical" 
dynamo could generate and amplify the regular magnetic field 
to the equilibrium field strength within ^1.5 Gyr (at z ^ 3.5; 
Fig. 121). At z ■X, 2, the UV background intensity decreases 
(Bianchi et al. 2001), resulting in the formation of a thin disk 
in a small galaxy that preserves the strength and ordering of the 
existing regular field. Although the ordering timescale of dwarf 
galaxies is longer (Fig. [3]), fully coherent regular magnetic 
fields were generated at earlier epochs (z ~ 1) because of the 
smaller sizes of dwarf galaxies. 

Due to the hierarchical clustering, it is expected that more 
massive protogalaxies formed at later epochs. In this scenario, 
the first disks should have formed in MW-type galaxies (see 
Fig. m and, hence, the time to reach magnetic-field strengths 
of ~ 10"'' G in MW-type and giant galaxies would be com- 
parable, while the ordering of the regular fields would have 
been completed earlier in MW-type galaxies. If dwarf galaxies 
formed earlier (z ^ 10), the regular field would have reached 
the equilibrium magnetic-field strength earlier and would have 
been ordered even earlier, at z ~ 1 (Fig.[3]l. 

6.4. Influence of star formation and mergers 

Star formation. The evolution of regular magnetic fields de- 
pends on the star-formation rate (SFR) and the disk parameters. 
Models of disk formation demonstrated that star formation is a 
fundamental parameter. It can be triggered in isolated galaxies 
by gravitational instability, in interacting galaxies by minor and 



major mergers, by tidal forces, leading to the compression of 
the gas, and by interactions of high velocity HI diffuse clouds 
(e.g. Kennicutt et al. 1987; Combes 2005). Star formation is 
more efficient at high redshifts because of the high merging rate 
and more gas at z - 1.3 (Lefevre et al. 2000; Ryan et al. 2008), 
and in more massive disk galaxies (Governato et al. 2007). 

Major mergers of gas-rich galaxies may trigger high SFR, 
causing high velocity turbulence of the ionized gas, which in 
turn can suppress the mean-field dynamo in the thick disk if the 
dynamo number is |Dd| < ~ 100 - 300 and in the thin disk 
if |Dc| < Dc ~ 7 (see Sect. 123] and IZTT). The latter restriction 
places an upper limit on the turbulent velocity of the gas: 

y ^ l.lQ/i= 11 kmsec"', (26) 

where Q = 20kms ' kpc ' and h - 0.5 kpc. The turbulent ve- 
locity (or characteristic velocity dispersion of the gas) and the 
SFR correlate positively for nearby galaxies (Dib et al. 2006). 
Using this relation and Eq. ( |26l ), one derives an upper limit of 
SFR $ 20 (in units of Mq yr"') up to which the action of the 
large-scale dynamo is possible. This indicates that the merging 
history of the galaxy had an important influence on the evo- 
lution of magnetic fields. In this scenario, regular fields were 
generated primarily at later epochs, z < 4 (see Fig.[T]i. 

Major mergers were rare but could alter the morphology 
of a spiral galaxy and destroy its regular magnetic field. After 
the merger, there were two possibilities, firstly, to form a disk 
galaxy with spheroidal or bulge component, or secondly, to 
form an elliptical galaxy without a disk. The SFR during the 
major merger 2x lO"^ yr) was ~ 50 times higher than in the 
case of an isolated galaxy (di Matteo et al. 2007; Bournaud 
et al. 2007). This high SFR produced an increase in total 
magnetic-field strength as shown by the peak of the equipar- 
tition magnetic-field strength at about 5 Gyr in Fig.[T] 

Disk galaxies, which survived after a gas-rich major 
merger, formed a thin-disk component immediately after the 
merger event and needed a; 1 .5 Gyr to amplify the regular mag- 
netic field to the equilibrium level by the mean-field dynamo 
(Fig. [TJ and ^ 8 Gyr to generate a coherent magnetic field on 
the length scale of the galaxy size (Fig.[3]i. 

Minor mergers were more frequent and may also have al- 
tered the morphology (spiral into elliptical, spiral to spheroid), 
increased the size and thickness of the disk, and controlled the 
star formation rate (density, turbulence velocity). Simulations 
of the formation of thick disks by minor mergers of satellites 
with a thin disk showed that mergers with 10-20% mass of the 
mass of the host galaxy resulted in the formation of thick disks. 
The disk was not destroyed but considerably heated up, tilted, 
and flared (Villalobos et al. 2008). The scale length increased 
slightly with respect to the original host disk, while the scale 
height became up to four times larger, depending on the incli- 
nation of the satellite's orbit. A larger disk height and radius 
increased both the dynamo timescale (Eq. [TJ and the ordering 
timescale (Eq. [T3T l, and shifted the formation of regular mag- 
netic fields to later epochs. 

Multiple minor mergers. Depending on the mass ratio of 
galaxies and the number of minor mergers, the disk could have 
been preserved, forming a spheroidal component, or destroyed. 
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forming multiple spheroids or an elliptical galaxy (Bournaud 
et al. 2007). The disk-formation model presented by Stringer 
& Benson (2007) illustrated that the resulting disk of an iso- 
lated galaxy and a galaxy formed by minor mergers have a very 
similar evolutionary behaviors in which the main governing pa- 
rameter is the infalling gas supplied from the hot halo. 

7. Discussion 

7.1. Observational tools 

The radio - far-infrared correlation provides a powerful tool 
for detecting star formation in distant galaxies with the help of 
ground-based radio observations. However, its validity depends 
critically on the magnetic-field strength coupled with the inter- 
stellar gas density and, hence the star-formation rate (Helou 
& Bicay 1993; Lisenfeld et al. 1996; Niklas & Beck 1997). 
In galaxies with recent starbursts, the radio continuum emis- 
sion can be deficient with respect to the far-infrared emission 
(Roussel et al. 2003). Either the cosmic rays or the magnetic 
field or both are generated with some time delay with respect 
to the dust-heating radiation field. 

The total magnetic field can be measured by the observed 
total power radio emission, corrected for the thermal fraction of 
a galaxy. The regular magnetic field can be traced by polarized 
synchrotron emission and by Faraday rotation. While Faraday 
rotation is an unambiguous signature of coherent regular fields, 
polarized emission can also emerge from "anisotropic" fields, 
which are the result of compressing or shearing of isotropic tur- 
bulent fields. This process has been observed in several galax- 
ies of the Virgo cluster (Vollmer et al. 2007; Wezgowiec et al. 
2007). 

7.2. Comparison with observations of galaxies 

Regular magnetic fields, as traced by polarized synchrotron 
emission and Faraday rotation, were detected in the disks of 
nearby spiral galaxies, with spiral patterns and amplitudes of 
up to 15 fiG, and were generally strongest in interarm regions 
(Krause 1990; Beck 2005). Mean-field dynamo theory, in spite 
of its simplifications, reproduces successfully the basic spiral 
pattern of magnetic fields observed in nearby spiral galaxies 
as well as in barred galaxies (Beck et al. 1996; Kulsrud 1999; 
Riidiger & HoUerbach 2004; Moss et al. 2007). The prominent 
example is M31, which hosts a coherent axisymmetric field 
within the emission "ring" at about 10 kpc radius (Fletcher et 
al. 2004). However, regular fields with a coherence scale simi- 
lar to the galaxy size are rare among nearby galaxies. 

In contrast to the turbulent dynamo, the mean-field dynamo 
needs far more time to generate a regular field. Fields with 
a coherence length of about 1 kpc can be expected after 1- 
2 Gyr, at z ^ 3 - 4 (Fig.O, but the generation time for regular 
fields with a coherence length of the galaxy size is several times 
longer and comparable to the galaxy age (Fig.|3]l. The ordering 
timescale for a 10 kpc coherence scale such as that in M 31 is 
about 10 Gyr (Fig.O. The coherent field of M 31 is an indica- 
tion that this galaxy has not experienced from a major merger 
during the past 10 Gyr. 



Most other nearby galaxies observed so far have more com- 
plicated field patterns than in M31 (Beck 2005). The long 
timescale required for field ordering (Fig.|3]l is consistent with 
these observations: most galaxies probably did not yet reach 
full field coherence, either due to their size or due to major 
merging events. The lifetime of galaxies larger than 15 kpc is 
insufficient to generate a fully coherent field. If a major merger 
in a MW-type galaxy occurred less than about 5 Gyr ago, the 
recovery of the fully coherent regular field is not yet completed. 

The results of this paper suggest a correlation between co- 
herence scale and galaxy size. The only direct indication of a 
regular field in a distant galaxy was derived from Faraday ro- 
tation in an intervening galaxy at z=0.395 measured against a 
background quasar (Kronberg et al. 1992). Since the pattern ex- 
hibits reversals and hence is incoherent on the galaxy scale, this 
result is consistent with our results. Indirect evidence of regu- 
lar fields in foreground spiral systems beyond z - I has been 
provided by Faraday rotation probes of distant radio sources 
(Kronberg et al. 2008; Bernet et al. 2008). They found that the 
scatter in the intrinsic RMs of spirals is ~ 20 - 50 rad irT^, 
which implies a coherence of the field on scales of about half 
the galaxy size (in the case of no field reversals). This coher- 
ence scale is consistent with the ordering scale of magnetic 
fields derived in this work (see Fig. |3]l for MW-type galaxies 
in the redshift range z ~ I - 2. 

In the LMC, another axisymmetric field was detected, 
which was found to have properties conflicting with dynamo 
theory (Gaensler et al. 2005). Since the LMC has a disk, al- 
though this rotates slowly, the mean-field dynamo can indeed 
operate and create a global coherent field on a scale of 4 kpc 
(radius of the LMC) within about 10 Gyr (Fig. [3). Hence, dwarf 
galaxies are prime candidates for hosting fully coherent regu- 
lar fields. On the other hand, dwarf galaxies are affected more 
significantly than larger galaxies by interactions. Simulations 
indicate that the LMC may have passed through the plane of 
the Milky Way more than 2 Gyr ago, or that it is on its way 
to a first passage through the MW (Kallivayalil et al. (2006), 
Piatek et al. (2008), and references therein). Since the regular 
field would be severely distorted during each passage, the ob- 
servation of a regular field supports the idea of no passages in 
the past. 

M 82 is the prototypical and most nearby starburst galaxy, 
excited by a tidal interaction with M 8 1 . Its total magnetic field 
is strong, about 50 yuG (Klein et al. 1988), and probably in 
equipartition with the other components of the wind outflow. 
The radio emission from the galaxy disk appears to be almost 
completely depolarized (Reuter et al. 1994). The radio halo 
is highly polarized but Faraday rotation is small, indicating 
anisotropic turbulent fields that are sheared and compressed 
in the rapidly moving gas outflow. This confirms our model 
of a strong and dominant turbulent field in starburst galaxies 
(Fig. [TJ. Observations of higher frequency, however, also de- 
tected polarized emission in the disk (Wielebinski 2006), indi- 
cating that some regular disk field may have survived the tidal 
interaction, which might not have been the case after a galaxy 
merger. To search for regular fields in the disks and halos of 
starburst galaxies, new observations of higher resolution are re- 
quired. 
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From Zeeman splitting in the HI absorption line against, 
a background quasar (Wolfe et al. 2008) a 84 ± 9 yuG field 
with a coherence length of > 200 pc was discovered in a low- 
metallicity galaxy at z = 0.692. A turbulent dynamo can gen- 
erate such a strong turbulent field in a starburst galaxy trig- 
gered by a major merger (see Fig. [U. However, this field can- 
not produce a Zeeman signal because of line-of-sight integra- 
tion through the whole galaxy, unless the observed field is lo- 
cated in a single dense gas cloud or filament. Alternatively, the 
observed field could originate in a shear or compression layer 
between interacting galaxies, although verification of this pos- 
siblity is beyond the scope of this paper. 

7.3. Detection of synchrotron emission from distant 
galaxies with SKA 

One major objective of SKA observations will be the detection 
of radio continuum emission from distant galaxies to study the 
star-formation history in the early Universe (van der Hulst et 
al. 2004). The tight radio - far-infrared correlation in galaxies 
implies that radio synchrotron emission is an excellent tracer 
of star formation in galaxies (Condon et al. 1991), at least to 
distances of z ^ 3 (Seymour et al. 2008). However, its applica- 
tion to even higher redshifts depends crucially on the existence 
of magnetic fields at the equipartition level with turbulent gas 
motions. 

The timescale for field generation in young galaxies con- 
strains the largest distance to which SKA can detect star- 
forming galaxies in synchrotron emission. We have shown in 
this paper that, even in the case of very weak seed fields of 
10 G, the turbulent small-scale dynamo, with the help of 
virial turbulence, can amplify turbulent fields to the level of 
equipartition with turbulent energy density within ^ 3 x 10^ 
years; strong fields should therefore exist in all star-forming 
galaxies at z ^ 10 (Fig. [1]) and the radio - far-infrared cor- 
relation should be valid for z ^ 10. Hence, radio continuum 
emission can be used as a tracer of star formation in the early 
Universe. 

Star-formation rate and total magnetic-field strength are re- 
lated nonlinearly (Niklas & Beck 1997; Chyzy et al. 2007; 
Krause 2008), so that young starburst galaxies should have 
much stronger fields and a higher radio - far-infrared ratio than 
normal galaxies. The SKA and its pathfinder telescopes will 
investigate this relation in more detail. 

The SKA will detect radio continuum emission from spi- 
ral galaxies similar to the MW to z - 3 and to even larger 
distances for more massive and starburst galaxies. The predic- 
tion from dynamo models presented in this paper can be tested 
with the help of polarization and Faraday rotation data from the 
SKA radio continuum and RM surveys (Gaensler et al. 2004). 
Systematic observations of galaxies at z = 1-5 are expected to 
envisage the early stages of galactic magnetic-field evolution. 

The first regular magnetic fields of contemporary disk 
galaxies were formed at z ~ 4, when the h/R ratio was 
higher than today. Therefore, their field patterns may have been 
substantially different from those of local galaxies (Moss & 
SokoloflF 2008), and their coherence scales were still smaller 



than their sizes. Our model predicts that dwarf and giant galax- 
ies have generate regular fields of several fiG strength by z ~ 4, 
while MW-type galaxies should not host these fields before 
z ~ 3. An important implication is that polarized radio disks 
and Faraday rotation are expected to exist in all galaxies at high 
redshifts (z ^ 3). Even if these galaxies cannot be resolved by 
the SKA beam, the integrated radio emission should be signif- 
icantly polarized for moderate and high inclinations (Stil et al. 
2008). A statistical analysis of a significant sample of galaxies 
observed by the SKA surveys would be able to test the valid- 
ity of the predictions of this paper. The detailed field pattern 
can be measured with the help of a dense RM grid of polarized 
background sources (Stepanov et al. 2008). 

Major merging events destroy the field regularity, while 
the total equipartition field strength is increased over a certain 
period (see Fig. [T]and Sect. I6.4l l. Furthermore, violent forces 
during the merging process can compress the turbulent field 
and lead to strong polarized emission. Hence, the polarized in- 
tensity increases, although the large-scale coherent field will 
be weaker than before. Faraday rotation measurements are re- 
quired to distinguish coherent from anisotropic fields. 

Some fraction of galaxies will be observed in a post-merger 
or starburst phase when the regular field has been temporar- 
ily destroyed and has not yet recovered (Fig. [U. A correlation 
between Faraday rotation and signatures of a recent merger is 
expected and should be investigated, whereas the total and po- 
larized emission should be higher due to enhanced turbulence 
and a fast outflow. Our estimates constrain the time since the 
peak of the starburst and will be useful for observing distant, 
starburst galaxies. 

Outflows from starburst galaxies contribute to the magne- 
tization of the intergalactic medium (IGM) (Kronberg et al. 
1999). Estimating the importance of this effect will require re- 
liable statistical data about starburst frequency, e.g. from radio 
observations with SKA. 

8. Conclusions 

Studying the evolution in magnetic fields of galaxies is impor- 
tant for interpreting future radio synchrotron observations with 
the planned Square Kilometre Array (SKA). We have used the 
dynamo theory to derive the timescales of amplification and or- 
dering of magnetic fields in disk and quasi-spherical galaxies. 
This has provided a useful tool in developing a simple evolu- 
tionary model of regular magnetic fields, coupled with mod- 
els describing the formation and evolution of galaxies. In the 
epoch of dark matter halo formation, seed magnetic fields of 
~ 10"'** G strength were generated in protogalaxies by the 
Biermann battery. Turbulence in the protogalactic halo gener- 
ated by thermal virialization could have driven the turbulent 
(small-scale) dynamo and amplify the seed field to the equipar- 
tition level of a; 20 jjG within a few 10** yr. In the epoch of disk 
formation, the turbulent field served as a seed for the mean-field 
(large-scale) dynamo developed in the disk. 

- We defined three characteristic timescales for the evolution 
of galactic magnetic fields: one for the amplification of the 
seed field, a second for the amplification of the large-scale 
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regular field, and a third for the field ordering on the galac- 
tic scale. 

- Galaxies similar to the Milky Way formed their disk atz~ 
10. Regular fields of equipartition (several juG) strength and 
a few kpc coherence length were generated within 2 Gyr 
(until z ~ 3), but field ordering up to the coherence scale of 
the galaxy size took another 6 Gyr (until z ~ 0.5). 

- Giant galaxies had already formed their disk at z ^ 10, al- 
lowing more efficient dynamo generation of equipartition 
regular fields (with a coherence length of about 1 kpc) un- 
til z ~ 4. However, the age of the Universe is too young 
for fully coherent fields to have akeady developed in giant 
galaxies larger than about 15 kpc. 

- Dwarf galaxies formed even earlier and should have hosted 
fully coherent fields atz~ 1 . 

- Major mergers excited starbursts with enhanced turbulence, 
which in turn amplified the turbulent field, whereas the reg- 
ular field was disrupted and required several Gyr to recover. 
Measurement of regular fields can serve as a clock for mea- 
suring the time since the last starburst event. 

- Starbursts due to major mergers enhance the turbulent field 
strength by a factor of a few and drive a fast wind outflow, 
which magnetizes the intergalactic medium. Observations 
of the radio emission from distant starburst galaxies can 
provide an estimate of the total magnetic-field strength in 
the IGM. 

This evolutionary scenario can be tested by measurements 
of polarized synchrotron emission and Faraday rotation with 
the SKA. We predict an anticoiTelation at fixed redshift be- 
tween galaxy size and the ratio between ordering scale and 
the galaxy size. Weak regular fields (small Faraday rotation) 
in galaxies at z ^ 3, possibly associated with strong anisotropic 
fields (strong polarized emission), would be signatures of major 
mergers. Undisturbed dwarf galaxies should host fully coherent 
fields, giving rise to strong Faraday rotation signals. 
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